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Abstract 



The definition and derivation of flavor neutrino states in the frame- 



> : 
o . 

work of standard Quantum Field Theory is reviewed, clarifying some 
^ . subtle points. It is shown that a flavor neutrino state that describes 

(Nj \ a neutrino produced or detected in a charged- current weak interaction 

process depends on the process under consideration and is appropriate 
for the description of neutrino oscillations as well as for the calculation 
of neutrino production or detection rates. Hence, we have a consistent 
framework for the description of neutrino oscillations and interactions 
in neutrino oscillation experiments. The standard flavor neutrino states 
are obtained as approximations which describe neutrinos in experiments 
that are not sensitive to the dependence of neutrino interactions on the 
neutrino mass differences. It is also shown that the oscillation proba- 
bility can be derived either through the usual light-ray time = distance 
Qh! approximation or through an average of the space-time dependent os- 

cillation probability over the unobserved propagation time. 
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1 Introduction 

The standard theory of neutrino oscillations was derived in 1975-76 [TJ [2], [3j 0] under the 
assumption that a flavor neutrino u a , produced or detected in a charged-current weak 
interaction process together with a charged lepton with flavor a (a = e, /x, r) is described 
by the standard flavor state 

k) = I>4h>, (i.i) 

k 

where U is the unitary mixing matrix of the neutrino fields and \v}.) are the Fock states 
of the massive neutrino fields, with corresponding masses (see the reviews in Refs. [HJ 

ElME!). 

In Ref. [10] it has been shown that the flavor state in Eq. (11 .11) is not a quantum of the 
flavor field u a (it is not annihilated by v a if the neutrino masses are taken into account). 
However, it was shown that the flavor state in Eq. (11.11) describes a flavor neutrino v a in 
the realistic ultrarelativistic approximation (m& <C E, where E is the neutrino energy) 

ma- 

Later [TT] . it has been discovered that it is possible to construct a Fock space of 
flavor states, which allows an alternative description of neutrino oscillations (see also 
Refs. [El UHl HH [151 H6l HZl [TBI [ini EQ]). However, this construction suffers from mys- 
terious features (the number of flavor Fock spaces is infinite, depending on arbitrary 
mass parameters; each flavor vacuum and the corresponding ladder operators are time- 
dependent; the ladder operators at different times do not satisfy the canonical anticom- 
mutation relations) and it has been shown that the flavor Fock states cannot be applied 
to the calculation of interaction processes [2TJ, [22] . 

Let us emphasize that the inapplicability of flavor Fock space theories to the cal- 
culation of neutrino interactions is a crucial shortcoming for their applicability to the 
description of neutrino oscillations, because any oscillation experiment involves the pro- 
duction and detection of neutrinos. 

It is then natural to ask if the standard theory of neutrino oscillations is appropriate 
for the description of neutrino interaction processes. In other words, are the standard 
flavor states in Eq. (11.11) appropriate for the calculation of the neutrino production and 
detection rates? 
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In this paper we review the derivation the flavor neutrino states in the framework 
of standard Quantum Field Theory, clarifying some subtle issues. We will show that 
the flavor neutrino states are appropriate for the description of neutrino production and 
detection, as well as for the description of neutrino oscillations^. We will see that the 
flavor neutrino states reduce to the standard flavor states in Eq. (11.11) in the case of 
experiments which are not sensitive to the difference of the neutrino masses. In this case, 
we recover the standard oscillation probability. 

In this paper we adopt the so-called "plane wave approximation" , in which the massive 
neutrino components of a flavor state are described by plane waves, as in the standard 
approach (see the reviews in Refs. [SI El El El Ej)- However, for the justification of some 
assumptions in the derivation of the oscillation probability, we will need to take into 
account the wave packet character of propagating massive neutrinos [25l EHJ E7J EHJ [29], 
EOl EH [321 E31 EH [351 [361 EZl EHl E9], albeit without need of a specific model. 

The plan of the paper is as follows. In section [2] we review the standard derivation of 
the neutrino oscillation probability, highlighting the underlying assumptions. In section El 
we present a derivation of the flavor neutrino states in the framework of the standard 
Quantum Field Theory. In section 0] we show that the flavor neutrino states are appro- 
priate for the description of neutrino production and detection. In section El using these 
flavor states, we derive the probability of neutrino oscillations, taking into account the 
general possibility that different massive neutrinos may have different momenta as well 
as different energies [101 123 HD H21 EH]- In section E] we show that the corrections to 
the standard oscillation phase due to violations of the light-ray time = distance approx- 
imation and the one-dimensional propagation approximation, adopted in section El are 
negligible. In section [7] we show that the standard oscillation probability can also be 
derived through the average over the unobserved propagation time of the space-time de- 
pendent oscillation probability, taking into account the general properties of the massive 
neutrino wave packets. 



2 Standard Derivation of the Neutrino Oscillation 
Probability 



where v a i,{x) are the left-handed flavor neutrino fields, v^[x) are the left-handed massive 
neutrino fields and U is the unitary mixing matrix (see the reviews in Refs. [S, El [71 El E]). 
Since a flavor neutrino v a is created by v aL {x) in a charged- current weak interaction 
process, in the standard plane- wave theory of neutrino oscillations PQ El El HI E]> ^ * s 
assumed that u a is described by the standard flavor state in Eq. ( 11. ip . which has the 
same mixing as the field 

1 In this review we consider neutrino oscillations on vacuum. We do not consider the more complicated 
case of neutrino oscillations in matter (23 [21] . 



Neutrino oscillations are a consequence of neutrino mixing: 





(2.1) 



A 
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The massive neutrino states \v k ) have definite mass m k and definite energy E k . Hence, 
they evolve in time as plane waves: 

i M*)> = J% \vk(t)) = Ek \u k {t)) =}► \u k (t)) = e~ iEkt \u k ) , (2.2) 

where is the free Hamiltonian operator and \v k it = 0)) = \v k ) (all the massive 
neutrinos start with the same arbitrary phase). The resulting time evolution of the flavor 
neutrino state Eq. (II. ip is given by 



K(t)) = J2K k e~' lEkt W k )= £ [^Kke-^U.AW,). (2.3) 




Hence, if the mixing matrix U is different from unity (i.e. if there is neutrino mixing), 
the state \u a (t)), which has pure flavor a at the initial time t = 0, evolves in time 
into a superposition of different flavors. The quantity in parentheses in Eq. (12.31) is the 
amplitude of v a — > up transitions at the time t after v a production. The probability of 
v a transitions at the time t = T of neutrino detection is given by 



2 

P Va ^JT) = \(vp\v a (T))[< " 



J2K k e-^ T U,k 

k 



-J^KkU^U^e-^-^. 

(2.4) 

One can see that P Va _ V/s (T) depends on the energy differences E k — Ej. In the standard 
theory of neutrino oscillations it is assumed that all massive neutrinos have the same 
momentum p. Since detectable neutrinos are ultrarelativisticd, we have 



I 7T)^ A 77)? 

E k = ^f + ml^E + ^ =► E k -E, = ^f, (2.5) 

where Am^ = ml — m| and E = \p\ is the energy of a massless neutrino (or, in other 
words, the neutrino energy in the massless approximation). In most neutrino oscilla- 
tion experiments the time T between production and detection is not measured, but 
the source-detector distance L is known. In this case, in order to apply the oscillation 
probability to the data analysis it is necessary to express t as a function of L. Consid- 
ering ultrarelativistic neutrinos, we have T ~ L, leading to the standard formula for the 
oscillation probability: 

P^(L, E) = Y1 Kk Upk U aj U; 3 exp f-i ■ (2.6) 

Summarizing, there are three main assumptions in the standard theory of neutrino 
oscillations: 

(Al) Neutrinos produced or detected in charged- current weak interaction processes are 
described by the flavor states in Eq. (11.11) . 



2 It is known that neutrino masses are smaller than about one eV (see the reviews in Refs. [151 144) ) . 
Since only neutrinos with energy larger than about 100 keV can be detected (see the discussion in 
Ref. [32j). in oscillation experiments neutrinos are always ultrarelativistic. 
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(A2) The massive neutrino states \i/k) in Eq. ( 11. ip have the same momentum (" equal- 
momentum assumption"). 

(A3) The propagation time is equal to the distance L traveled by the neutrino between 
production and detection ("time = distance assumption"). 



In the following we will show that the assumptions (Al) and (A3) correspond to ap- 
proximations which are appropriate in the analysis of current neutrino oscillation ex- 
periments. Instead, the equal-momentum assumption (A2) is not physically justified 



, j2],[33] , as one can easily understand from the application of energy- momentum 
conservation to the production process^]. However, in section we will show that the as- 



sumption (A2) is actually not necessary for the derivation of the oscillation probability 



if both the evolutions in space and in time of the neutrino state are taken into account. 



3 Flavor Neutrino States 

The state of a flavor neutrino v a is defined as the state which describes a neutrino pro- 
duced in a charged-current weak interaction process together with a charged lepton 
or from a charged lepton £~ = e ± ,fi ± ,r ± for a = e,/i,r, respectively), or the state 
which describes a neutrino detected in a charged-current weak interaction process with 
a charged lepton £~ in the final state. In fact, the neutrino flavor can only be mea- 
sured through the identification of the charged lepton associated with the neutrino in a 
charged-current weak interaction process. 

Let us first consider a neutrino produced in the generic decay process 

Pi->P F + £ + i/ a , (3.1) 

where Pi is the decaying particle and Pf represents any number of final particles. For 
example: in the pion decay process 

7T+ -> //+ + ^ , (3.2) 

we have Pj = n + , P F is absent and a = \i\ in a nuclear (3 + decay process N(A,Z) — > 
N(A, Z-l) + e + + u e we have Pi = N(A, Z), P F = N(A, Z - 1) and a = e. The following 
method can easily be modified in the case of a v a produced in the generic scattering 
process £~ + Pi — > Pf + v a by replacing the £+ in the final state with a £~ in the initial 
state. 

The final state resulting from the decay of the initial particle Pi is given by 

l/>=S|Pi>, (3.3) 

where S is the S- matrix operator. Since the final state |/) contains all the decay channels 
of Pi, it can be written as 

|/) = ^^I^C,Pf) + ... , (3.4) 

k 

3 A different opinion, in favor of the equal-momentum assumption, has been recently expressed in 
Ref. [IS]. On the other hand, other authors [1SJ U71 [35] advocated an equal-energy assumption, which 
we consider as unphysical as the equal-momentum assumption. 
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where we have singled out the decay channel in Eq. (13.11) and we have taken into account 
that the flavor neutrino v a is a coherent superposition of massive neutrinos v k . Since 
the states of the other decay channels represented by dots in Eq. ( 13. 4ft are orthogonal 
to \i/ k , £+,P F ) and the different states |z^fc, £+, P F ) are orthogonal and normalized, the 
coefficients A^ k are the amplitudes of production of the corresponding state in the decay 
channel in Eq. ( 13. lft : 

A v ak = P P |/) = P P |S|Pi>. (3.5) 

Projecting the final state in Eq. ( 13.41) over P F ) and normalizing, we obtain the flavor 
neutrino state [Hfl S21 El ED] 

k p >=fel^l 2 ) ' £^k>- ( 3 - 6 ) 

Therefore, a flavor neutrino state is a coherent superposition of massive neutrino states 
\v k ) and the coefficient A^ k of the A; th massive neutrino component is given by the ampli- 
tude of production of v k . Since, in general, the amplitudes A v ak depend on the production 
process, a flavor neutrino state depends on the production process. In the following, we 
will call a flavor neutrino state of the type in Eq. (13. 6p a "production flavor neutrino 
state" . 

Let us now consider the detection of a flavor neutrino v a through the generic charged- 
current weak interaction process 

i/„ + Di->D F + C, (3.7) 

where Di is the target particle and Dp represents one or more final particles. In general, 
since the incoming neutrino state in the detection process is a superposition of massive 
neutrino states, it may not have a definite flavor. Therefore, we must consider the generic 
process 

f + Di, (3.8) 

with a generic incoming neutrino state \v). In this case, the final state of the scattering 
process is given by 

|/) = S|i/,D I ), (3.9) 
This final state contains all the possible scattering channels: 

\f) = |D F , £-) + ..., (3.10) 

where we have singled out the scattering channel in Eq. (13.71) . We want to find the 
component 

|i/ a ,Di) = 5]^ fc |^,Di> (3.11) 

k 

of the initial state \u, Di) which corresponds to the flavor a, i.e. the component which 
generates only the scattering channel in Eq. ( 13.71) . This means that |D F , £~) = S \u a , Di). 
Using the unitarity of the mixing matrix, we obtain 

|i/ a ,Di)=St|D F ,0- (3-12) 
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From Eqs. ( 13.111) and (13. 12j) . the coefficients A® k are the complex conjugate of the am- 
plitude of detection of v k in the detection process in Eq. ( 13. TP : 

^ = (i/ Jfe ,D I |St|D F ,0- ( 3 - 13 ) 

Projecting \i/ a , Di) over |Di) and normalizing, we finally obtain the flavor neutrino state 
in the detection process in Eq. (13.71) : 

iO = fei^i 2 ) ' XXm- ( 3 - 14 ) 

In the following, we will call a flavor neutrino state of this type a " detection flavor neutrino 
state" . 

Although the expressions in Eqs. (I3.6P and (13.141) for the production and detection 
flavor neutrino states have the same structure, these states have different meanings. 
A production flavor neutrino state describes the neutrino created in a charged- current 
interaction process, which propagates out of a source. Hence, it describes the initial 
state of a propagating neutrino. A detection flavor neutrino state does not describe a 
propagating neutrino. It describes the component of the state of a propagating neutrino 
which can generate a charged lepton with appropriate flavor through a charged-current 
weak interaction with an appropriate target particle. In other words, the scalar product 

A a = (3.15) 

is the probability amplitude to find a u a by observing the scattering channel in Eq. (13. 71) 
with the scattering process in Eq. 



In order to understand the connection of the production and detection flavor neutrino 
states with the standard flavor neutrino states in Eq. ( 11.11) . it is useful to express the 
S'-matrix operator as 

S = 1 - % I d 4 x J^cc{x) , (3.16) 



where we have considered only the first order perturbative contribution of the effective 
low-energy charged-current weak interaction Hamiltonian 

jebc(x) = ^jl(x)j"(x) 1 (3.17) 

where Gp is the Fermi constant. The weak charged current j p (x) is given by 



f(x) 



^(x) 7 p (l-7 5 )4(x) + //(x) 
E E ^ V (1 - 7 5 ) ia(x) + h"{x) , (3.18) 

Q=e,|i,r k 



where h p (x) is the hadronic weak charged current. The production and detection ampli- 

De written as 

A p ak = U* ak Ml k , A» k = U* ak M» k , (3.19) 



tudes A^ k and A® k can be written as 
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with the interaction matrix elements 

M p ak = ~i%J d A x (u k , £+\ W k {x) Y (1 - 7 5 ) IJx) |0) J P ^ (x) , (3.20) 

M ^ =i %S dAx H Wk{x) Y (1 " 75) Ux) 10 j °" Df * (x) • (3 - 2i) 

Here J Pl ^ Pp (x) and J° i_ * Df (x) are, respectively, the matrix elements of the Pi — > Pf 
and Di — > Dp transitions. 

Using Eq. (13. 19ft . the production and detection flavor neutrino states can be written 

as 



A/f D 

= E / Uak k> • (3-23) 

These states have a structure which is similar to the standard flavor states in Eq. (11. II) . 
with the relative contribution of the massive neutrino v k proportional to U* k . The ad- 
ditional factors are due to the dependence of the production and detection processes on 
the neutrino masses. 

In experiments which are not sensitive to the dependence of M^ k or M^ k on the 
difference of the neutrino masses it is possible to approximate 

M p ak z> M p a , AC ~ Ml . (3.24) 

In this case, since 



Y,\U ak \ 2 = l, (3.25) 



k 

we obtain, up to an irrelevant phase, the standard flavor neutrino states in Eq. (11.11) . 
which do not depend on the production or detection process. Hence, the standard flavor 
neutrino states are approximations of the production and detection flavor neutrino states 
in experiments which are not sensitive to the dependence of the neutrino interaction rate 
on the difference of the neutrino masses. 

In the following section H] we will show that the correct expressions for the production 
and detection flavor neutrino states are important in order to be able to describe, in a 
consistent framework, neutrino oscillations and neutrino production and detection. Then, 
in the next section [3] we will derive the neutrino oscillation probability starting from the 
production and detection flavor neutrino states. We will show that, with the appropriate 
approximations, the oscillation probability reduces to the standard one in Eq. (12. 6ft . which 
is derived from the approximate flavor neutrino states in Eq. (II. ip . 

4 Production and Detection Rates 

In order to measure v a z/g oscillations, it is necessary to calculate the neutrino pro- 
duction rate T a (E) of v a in the source and the detection cross section (Jp{E) of v@. The 
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number of transition events as a function of the distance L traveled by the neutrino 
between production and detection and the neutrino energy E is given by 

N a/3 (L, E) oc T a (E) P Va -+ Vj3 (L, E) ap(E) , (4.1) 

with a constant of proportionality which depends on the size and composition of the 
source and detector and on the running time of the experiment. From the measurement 
of N a p(L,E) and the knowledge of T a (E) and ap(E), the experimentalist infers the 
value of Pv a ^ l//3 (L, E), which gives information on the mixing parameters (elements of 
the mixing matrix and squared-mass differences) through Eq. (12. 6j) . 

Decay rates and cross sections are given by the incoherent sum over the different 
channels corresponding to different massive neutrinos [SH E21 ESI EH ES]- The reason 
is that massive neutrinos are the physical particles which propagate in space-time with 
definite kinematical properties. However, in section El we derived the oscillation probabil- 
ity P Va ^ V(3 (L,E) starting from the description of neutrinos through flavor states, which 
are coherent superpositions of massive neutrino states. Then it is natural to ask if the 
derivation of the oscillation probability is consistent with the calculation of decay rates 
and cross sections. In the following, we show that the description of neutrinos through fla- 
vor states leads to the correct expression for the decay rates and cross sections. Therefore, 
all quantities in Eq. (14.11) can be calculated in a consistent way. 

In order to be definite, we consider the general decay process in Eq. (13.11) . in which 
a flavor neutrino v a is produced. Any other process of neutrino production or detection 
can be treated in an analogous way. 

Using the flavor state in Eq. (13.61) and taking into account Eq. (13.51) . the amplitude 
of the general decay process in Eq. (13.11) is given by [50] 

K = <i£,£,Pf|S|Pi> = (eki 2 ) E^^' £ - p Fl s l p i) = ^Ei^i 2 - 

(4.2) 

Therefore, the decay probability is correctly given by an incoherent sum of the probabil- 
ities of production of different massive neutrinos, 

K| 2 = £l^| 2 - (4-3) 

i 

In other words, the coherent character of the flavor state in Eq. (13.61) is irrelevant for 
the decay probability, which can be obtained either using the flavor neutrino state in 
Eq. (13.61) or an incoherent mixture of massive neutrino states. The decay rate is then 
obtained by integrating each massive neutrino contribution to the decay probability over 
its phase space. 

Using the expression in Eq. (13. 19j) for the amplitude A^ k , the decay probability in 
Eq. (14.31) can be written as 

Kl 2 = £|tU 2 |.< fc | 2 , ( 4 - 4 ) 

k 

which is an incoherent sum of the probabilities of production of the different massive 
neutrinos weighted by \U ak \ 2 [5T | 1521 [53| [541 [55] . 
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Therefore, the flavor neutrino state in Eq. (13.61) leads to the correct decay rate for the 
general decay process in Eq. ( 13. ip . It is clear that this proof can easily be generalized to 
any charged- current weak interaction process in which flavor neutrinos are produced or 
detected. 

If an experiment is not sensitive to the dependence of Ai^ k on the different neutrino 
masses, it is possible to use the approximation in Eq. ( 13.24ft . In this case, using Eq. ( 13.251) . 
we obtain 

\Alf = \Ml\\ (4.5) 

If the scale of neutrino masses is negligible in comparison with the experimental resolution, 
the dependence of on the neutrino masses is negligible and the decay probability in 
Eq. (14.51) reduces to the standard decay probability for massless neutrinos. 

The decay probability in Eq. (14. 5 p can also be obtained starting from the standard 
flavor states in Eq. (11. ip . which are obtained from Eq. (13.61) through the approximation 
in Eq. ( 13.24j> . Indeed, in this case the decay amplitude is given by [56] 

Al = (u a , C, P F | S |Pi> = u <* A lk = E \Uak\ 2 Ml = Ml. (4.6) 

k k 

Let us remark, however, that in the case of an experiment which is sensitive to the 
dependence of M^ k on the different neutrino masses, a derivation of the decay amplitude 
starting from the standard flavor states would lead to an incorrect result. This is due to 
the fact that in this case the approximation in Eq. ( 13.241) is not valid and one must take 
into account the dependence of M^ k on the different neutrino masses in the definition of 
the flavor states. 



5 Neutrino Oscillations 

Let us consider a neutrino oscillation experiment in which u a —>■ u@ transitions are studied 
with a production process of the type in Eq. (13.11) and a detection process of the type in 
Eq. ( 13. 7p . In this case, the produced flavor neutrino u a is described by the production 
flavor state \u*) in Eq. ( 13. 6p . If the neutrino production and detection processes are sepa- 
rated by a space-time interval (L,T), the neutrino propagates freely between production 
and detection, evolving into the state 

HL,T)) = e-^ T ^ 1 } (5.1) 

where P° and P are, respectively, the energy and momentum operators. This is the 
incoming neutrino state in the detection process. The amplitude of the measurable u a — ► 
up transitions is given by the scalar product in Eq. ( 13. 151) : 

A„ a ^(L,T) = (up\v(L,T)) = (^\e-* 0T+ * rL \vl) , (5.2) 

with the detection flavor state \u9) in Eq. (13.141) . 

Since the massive neutrinos have definite kinematical properties (energy and momen- 
tum), we have, in the plane wave approximation, 

P"k> =p£k>> (5-3) 
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with 

p° k = E k = ^\p k \ 2 +m 2 k . (5.4) 
Using the normalization (vk\vj) = Skj, we obtain the flavor transition amplitude 

^,&T) = foj&A fe|./£| a ) ' ^^e"^^, (5.5) 

Notice that the consideration of the space-time interval between neutrino production and 
detection allows one to take into account both the differences in energy and momentum 
of massive neutrinos [HI [23 HH H21 E2! • 

Let us consider the simplest case in which all massive neutrino momenta p k are aligned 
along L. This assumption is reasonable, because all massive neutrino components are 
created in the same microscopic production process and detected in the same microscopic 
detection process, after propagation through the large macroscopic space interval L. In 
section [5] we will show that possible deviations from this assumption do not lead to any 
observable effect. 

In this "one-dimensional approximation", the transition amplitude depends on L = 
\L\ and T: 

where p k = \p k \. 

In oscillation experiments in which the neutrino propagation time T is not measured, 
it is possible to adopt the light-ray T = L approximation, since neutrinos are ultrarel- 
ativistic (the effects of possible deviations from T = L are shown to be negligible in 
section [6]). In this case, the phase in Eq. (15.61) becomes 

- E k T + Pk L = - (E k - Pk )L = — ^ L = - k L ~ % L , (5.7) 

1 ; E k + Pk E k +p k 2E 1 ; 

where E is the neutrino energy neglecting mass contributions. Equation ( 15. 7ft shows that 
the phases of massive neutrinos relevant for the oscillations are independent from the 
values of the energies and momenta of different massive neutrinos [101 EZl HD 1121 
because of the relativistic dispersion relation in Eq. (15. 4p . In particular, Eq. (15.71) shows 



that the equal-momentum assumption (A2) in section El adopted in the standard deriva- 



tion of the neutrino oscillation probability, is not necessary in an improved derivation 
which takes into account both the evolutions in space and in time of the neutrino state. 
The probability of v a —>■ Pp transitions in space is given by 



u a ^(L,E) = l^l 2 ) (E l^gl 2 ) ^Al k A^A**Al exp (- 



P., -„AL.E) = I > IAV.fl I > \A L l\- I > Al',.A L u. A v ''A u t; <>xi) I -i^ m W L 



2E 
(5. 
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Using the decomposition in Eq. (13.19j) . the oscillation probability in Eq. ( 15. 8ft 
be written as 



can now 



/ 



P Ua ^(L,E)= 

k,j 



Ml k Ml* 

Ei^riAcr 



\ 



/ 



\ 



^2\Upi\ 2 \M 



D \2 



x Kk Up k U aj Ul; exp 



' 2E 



(5.9) 



This probability has the same structure as the standard oscillation probability in Eq. ( 12.61) . 
with additional factors that take into account the effect of the neutrino masses in the 
production and detection processes. It is clear from Eq. (15.91) that these effects have an 
influence on the amplitude of the oscillations, but not on the phase, which coincides with 
the standard one in Eq. (12.61) . 

Since neutrinos in oscillation experiments are ultra-relativistic and the experiments 
are not sensitive to the dependence of neutrino interactions on the neutrino masses, the 
dependence of M^ k and M^ k on the neutrino masses can be neglected, leading to the 
approximation in Eq. (13.241) . In this case, the transition probability in Eq. (I5.9P reduces 
to the standard one in Eq. (12.61) . which can be obtained starting from the standard flavor 
states in Eq. ( 11.11) . As shown in section [31 the standard flavor states are obtained from the 
production and detection flavor states under the same approximations. Therefore, the 
standard flavor states are appropriate for the description of neutrino oscillation experi- 
ments in the plane wave approximation, as long as the dependence of the production and 
detection probabilities on the neutrino masses is negligible. These considerations justify 



the assumption (Al) in section El adopted in the standard derivation of the neutrino 



oscillation probability. 



6 Universality of the Oscillation Phases 

In the previous section [5l we derived the probability of neutrino oscillations under the 
assumptions that T = L and all massive neutrino momenta p k are aligned along L. In 
this section we will show that possible deviations from these assumptions do not affect 
in a significant way the oscillation phases measured in neutrino oscillation experiments, 
which are correctly given by the standard expression in Eq. ( 12.61) . 

The momentum of a massive neutrino created in a production process depends on the 
characteristics of the interaction, on the nature of the other particles taking part in the 
process and on the neutrino mass. Let us call p and E — \p\, respectively, the momentum 
and energy of a massless neutrino. The first order contribution of the neutrino mass to p k 
and E k must be proportional to ml, because the energy- momentum dispersion relation 
in Eq. ( 15.41) depends on m\. Therefore, in general, the momentum p k can be written to 
first order in m 2 k as 

2 

A-P-Cj^, (6-1) 



where 

f d Pk 



2E dm\ 



(6.2) 

m k =0 



12 



The value of the vector £ depends on the production process. However, as we will see 
in the following, the measurable oscillation phases are independent of £. Therefore, they 
are universal, i.e. independent of the specific nature of the neutrino production process, 
as well as detection. 

In the approximation in Eq. (jfi.ip . the dispersion relation in Eq. (15.41) implies that the 
energy of v k is given by 

Note that Eq. (16. ip implies that the directions of propagation of different massive 
neutrinos are slightly different if £ is not collinear with p. How this is possible can be 
illustrated in the case of the pion decay production process in Eq. (13.21) as follows. In 
the rest-frame of the pion, the energy of v k can be calculated from energy-momentum 
conservation to be given by 

m w / mf.\ mi 



where and m M are the masses of the pion and muon, respectively. Since, in the rest- 
frame of the pion, the neutrino and the muon are emitted back-to back, all the massive 
neutrinos are emitted in the same direction. Thus, £ is collinear with p, leading to 

2 2 
777 77? 

E k = E + (l-0^, Pk = \p k \^E-t^, (6.5) 

with 

777 / 771^ \ -» 1 / TV? \ 

Note that the different values of the momenta and energies of different massive neutrinos 
imply different corresponding values for the momentum and energy of the outgoing muon. 
Two different massive neutrinos can be produced coherently in the same decay process 
only if the outgoing muon has energy and momentum uncertainties which are larger, 
respectively, of the difference of the energies and momenta of the two massive neutrinos. 
These uncertainties must come from corresponding uncertainties for the pion. Hence, in 
a rigorous treatment, the pion and the muon must be described by wave packets [HH [35l 
I3T1 1361 13TI 1381 [39j [32]. This implies that also the massive neutrinos must be described 
by wave packets (2H1 [26l EH EHl EQl EH E21 E3] . The spatial extension of the massive 
neutrino wave packets explains how it is possible that the different massive neutrinos 
can be detected in the same interaction process even if their space-time trajectories are 
different. 

Let us now consider a boosted reference frame. If the frame is boosted in the direction 
of the emitted massive neutrinos, obviously they are still collinear. Instead, if the frame 
is boosted in another direction, collinearity is lost. For example, let us consider a frame 
boosted with velocity V in a direction which is orthogonal to the neutrino direction in 
the pion rest frame. In this frame, we have, in an obvious notation, 

P'k\\ = Pk\\ = Pk , p' k± = --j=== E k . (6.7) 
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Hence, the angle of propagation of Vk in the boosted frame, with respect to its direction 
in the rest frame, is given by 

tan^ = ^--^=fl + ^|l • (6.8) 



Since this angle depends on the neutrino mass, in the boosted frame different massive 
neutrinos are seen to propagate in slightly different directions. In fact, in this case, p 
and £ are not collinear, since 

P\\=P\\=E, p' ± = --j=L=E, tf, = $|=e, el=^=(l-0- (6-9) 

Since many experiments measure the oscillations of neutrinos produced by pion decay in 
flight (e.g. atmospheric and accelerator neutrino oscillation experiments), in which the 
neutrino direction is not constrained to be collinear with the pion direction, in order to 
be realistic, in this section we consider the general case of not-collinear p and £. 
Let us now calculate the phases 

<Pk = -E k T + p k -L, (6.10) 

in the transition amplitude in Eq. ( 15. 5ft considering a deviation of the time T between 
neutrino production and detection from the light-ray approximation T = L: 



T = W\ {1 + £Th (6 ' n) 



with the average velocity 



N 

Pk P 



N ^ E k E 

k=i K 




E 



m 2 
2E 2 



(6.12) 



where N is the number of massive neutrinos and 

N 

X 



' TrT, m l ( 6 - 13 ) 



m 2 



k=l 



is the average of the squared neutrino masses. The absolute value of the average velocity 
is given by the usual ultra-relativistic expression 

TTL 2 

\{v)\ ( 6 - 14 ) 

In Eq. (16. lip , we consider £y< 1, because neutrinos are ultra-relativistic and the devia- 
tion from T = L/\(v)\ cannot be larger than the size of the neutrino wave packets, which 
must be much smaller than L in order to observe the oscillations [2SI EZj- To first order 
in the small ratio m 2 /E 2 , Eq. (16.111) becomes 

L(l + e T )[l + ^- 2 \ ■ (6.15) 
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Now, we consider also deviations from the assumption that the momenta of all mas- 
sive neutrinos are aligned along L. As discussed above, the expression for the neutrino 
momenta in Eq. ( 16.11) already implies that, in general, the momenta of different massive 
neutrinos are not collinear. Moreover, we can consider a deviation from the collinearity 
of L and the average momentum 

1 N ~ 2 

k=l 



which we can write as 



^ + (6.17) 



with \el\ <C 1, for the same reasons of et <C 1- At zeroth order in the small ratio m 2 /E 2 , 
we have 



which implies, from Eq. 




(6.18) 



(6.19) 



From Eqs. (16. lj) . (16.15!) and (16. 19ft . for the difference Acpkj = 4>k — <fij °f the phases in 
Eq. (16.1 Op we obtain, at first order in the small quantities Et and \el\, 

AmhL 

A0 fcj * g|- + e kj , (6.20) 

with the contribution 



- < I, L ^ 



2£ 



(6.21) 



in addition to the standard oscillation phase in Eq. (I2.6p . However, since £y < 1 and 
\el\ 1, the contribution is non-negligible only for Avn^-LjlE ^> 1. But in this 
case oscillations are not measurable, since they are washed out by the average over the 
energy resolution of the detector (see Ref. [7]). In the case of AmljL/2E ~ 1, in which 
oscillations are measurable, Ekj is extremely small and can be safely neglected (e tSkj ~ 1), 
leading to the validity of the standard expression in Eq. (12. 6 p for the oscillation phases. 

Note the irrelevance for A0^j of the lack of collinearity of the trajectories of different 
massive neutrinos if Et = |?zj =0. It is due to the fact that the deviation from collinearity 
of Vk and Vj is proportional to Am 2 ,-/2E 2 (see the example in Eq. (16. 8ft ). Thus, it induces 
in the phase difference A(pkj effects of higher order in Am 2 k ^/2E 2 , which are completely 
negligible. 

In conclusion, in this section we have shown that possible small deviations from the 
light-ray approximation T = L and from the collinearity of the massive neutrino mo- 
menta p k and L are irrelevant for the measurable oscillation phases, which are universally 
independent from the specific characteristics of the process of neutrino production, as 



well as detection. Hence, the time = distance assumption (A3) in section [2} adopted in 



the standard derivation of the neutrino oscillation probability, is correct. 
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7 Time Average 



In section [5] we derived the probability of flavor neutrino oscillations as a function of the 
source-detector distance L through the light-ray T = L approximation. In the previous 
section [U] we have shown that the measurable oscillation phases are stable against possible 
small deviations from the light-ray approximation. Another way to obtain the oscilla- 
tion probability in space is through an average of the space-time dependent oscillation 
probability over the unobserved propagation time T [27], 

P^u,{L) oc j dTP„ a ^(L,T) , (7.1) 

where the constant of proportionality must be chosen in order to satisfy the conservation 
of probability constraint 

^P^(L) = J2 P ^u d (L,T) = 1. (7.2) 

P P 

The time average is done on the oscillation probability, not on the amplitude, since it is 
an average over the propagation times of different neutrinos in a beam, which contribute 
incoherently to the oscillation probability. 

Let us work, for simplicity, in the approximation in Eq. (13.241) for the production and 
detection matrix elements. From the one- dimensional approximation of the oscillation 
amplitude in Eq. (15. 6ft . we obtain the space-time dependent oscillation probability 

p^(l, t) = Kk Upk u aj u; 3 e -i(E k -E 3) T +l{Pk - P] )L {73) 

k,j 

The integration in Eq. (17.11) leads to 

P Va ^ (L) oc Kk Upk U aj W P3 e^~^ L 5(E k - E 3 ) . (7.4) 

k,j 

Hence, the time average implies that E k = Ej [39, 38J. However, from Eq. (I6.3P it is clear 
that in the plane wave approximation this equality is, in general, not possible. In order 
to allow such an equality, it is necessary to take into account the wave packet character 
of massive neutrinos, which implies an uncertainty in energy. 

In an accurate wave packet treatment, the massive neutrinos are described by super- 
positions of plane waves with different momenta [26j [271 [23 [291 EQJ SB E21 S3] . Since 
we are not interested in the wave packet effects, we adopt, for each massive neutrino z/^, 
an approximate wave packet with a distribution in energy ^(E^ — E k ), where Ej. is the 
average energy. In this case, the space-time dependent oscillation probability in Eq. (17. 3ft 
becomes 

P Ua ^ U0 (L, T) = U* k Up k U aj Upj 

x / dE k [ &E 3 <P(E k - E k ) <\>\E 3 - Ej) e -^ k -E 3 )T +l{Pk - P] )L ? (7 5) 
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with pk = \JE\ — mf. The integration over T in Eq. (17. ID leads to 

P Va ^{L) ^YsKkUpkU^U^ J dE <p(E — E k ) (jf(E — Ej) e i{pk ~ p ^ L , (7.6) 



with pk = tJE 2 — m\. In the ultrarelativistic approximation pk ~ E — m\/2E, we obtain 
P Va ^(L) ( xJ2 U * a kUpkU aj U; j J dE <j){E — Ek) (j>*(E — Ej) exp f-i^^j . (7.7) 

kj 

Now, we assume that the energy width of the wave packets is much larger than the 
difference between Ef. and Ej, 

4>(E - E k ) ~ <j>(E - Ej) ~ 0(£ - £) , (7.8) 

where E may be chosen as the average energy of the wave packet of a massless neutrino 
(or the average among the Ek's). This is the condition which allows the equality Ek = Ej 
without suppressing the oscillation probability. Since the energy width of the wave packets 
is inversely proportional to their spatial width and the difference \Ek — Ej\ is inversely 
proportional to the oscillation lengtb@, the condition in Eq. (17.81) is verified by spatial 
wave packets which are much smaller than the oscillation length. This is a necessary 
condition for the observation of neutrino oscillation which is satisfied in all experiments. 
Using the approximation in Eq. ( 17.81) . the oscillation probability in Eq. ( 17.71) becomes 

/ Am? L\ 

dE\cf>(E-E)\ 2 expl-i^-) . (7.9) 

Finally, we consider a sharply peaked wave packet, i.e. a wave packet <fi(E — E) 
with energy uncertainty much smaller than the average energy E. This is a realistic 
assumption, since in practice the energy uncertainties of the production and detection 
processes are much smaller than the neutrino energy. In this case, the phase Am\-L/2E 
is practically constant over the size of the wave packet when the distance L is of the order 
of the average oscillation length 4irE/ Amj,j. As we have already remarked in sectional 

this is a necessary condition for the observation of oscillations, since for L ^> AnE/Amlj 
the oscillations are washed out by the average over the energy resolution of the detector. 
Hence, using the normalization J dE \<j){E — E)\ 2 = 1, we finally obtain 

Pv a ^ (L) = Kk Upk U aj U* pj exp f-i • (7.10) 

kj V 2E J 

This is the standard oscillation probability in Eq. ( 12. 6ft . with E = E. 

Summarizing, in this section we have shown that the average of the space-time depen- 
dent oscillation probability over the unobserved propagation time T leads to the standard 
oscillation probability in space, taking into account the wave packet character of massive 
neutrinos, without need of a detailed wave packet model. 



4 Writing E k ~ E + (1 - £,)m 2 k /2E, in analogy with Eq. (|53]l . we have \E k - Ej\ oc \Am 2 kj \/2E 
1/L^ C , where LS C = 4nE/\Am%.j\ is the oscillation length. 
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8 Conclusions 



In this review we have presented a consistent framework for the description of neutrino 
oscillations and interactions, which is appropriate for the theoretical interpretation of 
neutrino oscillation experiments, in which neutrinos are produced by a source and de- 
tected after propagation over a macroscopic distance. We have shown that the flavor 
neutrino state that describes a neutrino produced or detected in a charged-current weak 
interaction process depends on the process under consideration and is appropriate for the 
description of neutrino oscillations as well as for the calculation of neutrino production 
and detection rates. A flavor neutrino state can be approximated with the standard ex- 
pression in Eq. ( 11. ip only for experiments which are not sensitive to the dependence of 
neutrino interactions on the different neutrino masses. This is the case of all neutrino 
oscillation experiments. 

We have also reviewed and clarified some subtle points concerning the derivation of the 
oscillation probability. In particular, we have shown that the oscillation probability can be 
derived through the usual light-ray time = distance approximation, because possible small 
deviations have negligible effects on the measurable oscillation phase. We have also shown 
that the oscillation probability can be derived in an alternative way through an average of 
the space-time dependent oscillation probability over the unobserved propagation time, 
taking into account the wave packet character of massive neutrinos, without need of a 
specific model. 
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